Dissipative Optomechanics in a Michelson— Sagnac Interferometer 
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Dissipative optomechanics studies the coupling of the motion of an optical element to the decay 
rate of a cavity. We propose and theoretically explore a realization of this system in the optical 
domain, using a combined Michelson-Sagnac interferometer, which enables a strong and tunable 
dissipative coupling. Quantum interference in such a setup results in the suppression of the lower 
motional sideband, leading to strongly enhanced cooling in the non-sideband-resolved regime. With 
state-of-the-art parameters, ground-state cooling and low-power quantum-limited position transduc- 
tion are both possible. The possibility of a strong and tunable dissipative coupling opens up a new 
route towards observation of fundamental optomechanical effects such as ponderomotive squeezing 
or nonlinear dynamics. Beyond optomechanics, the method suggested here can be readily transferred 
to other setups involving such systems as nonlinear media, atomic ensembles, or single atoms. 

PACS numbers: 42.50.Wk, 42.79.Gn, 07.10. Cm, 07.60.Ly 



Recent progress in the engineering of high-quality mi- 
cromechanical oscillators coupled to high-finesse cavity 
modes has paved the way towards sensing and con- 
trol of mechanical motion at the quantum limit [1-5]. 
The rapid developments in the field of optomechanics 
bear important implications for both applied and ba- 
sic science, ranging from applications in high-sensitivity 
metrology [6-8] and quantum information processing [9- 
13] to fundamental tests of quantum mechanics at large 
mass- and length-scales [14-16]. 

In the conventional paradigm of optomechanics the in- 
teraction of the mechanical oscillator with a cavity mode 
is dispersive in the sense that the cavity resonance fre- 
quency experiences a shift depending on the displacement 
of the mechanical oscillator arising from conservative ra- 
diation pressure or optical gradient forces. This coherent 
dispersive interaction has been employed for sideband- 
cooling to the quantum mechanical ground state [17, 18], 
as well as for the observation of optomechanical normal- 
mode splitting [19-21] and optomcchanically-induced 
transparency [22, 23]. The complementary paradigm of 
dissipative coupling — where the width Kc of the cavity 
resonance, rather than its frequency Wc, is dependent on 
the mechanical displacement x — was introduced very re- 
cently in a theoretical study [24] in the context of elec- 
tromechanics. This situation, rather unusual for cavity 
quantum electrodynamics, was shown [24] to give rise to 
remarkable quantum noise interference effects which dra- 
matically relax requirements for cooling to the ground 
state without sideband resolution; it also allows reaching 
the standard quantum limit (SQL) for the imprecision in 
position measurements. While it is already clear that a 
strong and tunable dissipative optomechanical coupling 
would greatly enrich the toolbox of optomechanics, its 
full significance for the quantum control of optomechani- 
cal systems, e.g., for ponderomotive squeezing [25], non- 
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FIG. 1. (a) Topology of the Michelson-Sagnac interferometer, 
(b) Effective cavity; the properties of M depend on M. 



linear dymanics [26, 27], or pulsed protocols [28], is yet 
to be explored. 

Unfortunately, an optomechanical setup having a strong 
dissipative coupUng in the absence of dispersive coupling 
has not been found. Consider, e.g., a Fabry-Perot in- 
terferometer (FPI) of length L and resonance frequency 
LUc = Tmc/L (n G W), with one movable ideal end-mirror 
and an input coupler of transmissivity r, such that the 
cavity linewidth is Kc = c]t]^/(4L). When the mirror 
moves, the cavity length changes, and with it both lOc 
and Kc. The corresponding shifts per zero-point fluctu- 
ation xq of the mirror oscillator, g^^ — {dujc/dx)xo and 
g^ ~ {dKc/dx)xo, quantify the strength of dispersive and 
dissipative coupUngs, respectively. Their ratio is thus 
given by the cavity's quality factor, goj/gn — Wc/kc ^ 1, 
such that the dispersive coupling dominates by far. A 
similar conclusion applies for a movable membrane cou- 
pled to a single mode of a FPL One may obtain 5^ ~ g^^ 
by coupling the membrane to multiple [29, 30] transverse 
modes, such that both types of coupling contribute to the 
dynamics; however in such a setup, one cannot 'switch 
off' the dispersive interaction to take advantage of the 
quantum noise interference effects present in a purely dis- 
sipative coupling. 

Here we show that in a Michelson-Sagnac interferon!- 



eter (MSI) with a movable membrane [31, 32], cf. Fig. 1, 
a strong and tunable optomechanical coupling can be 
achieved for which g^j — but where the effective dis- 
sipative coupling strength [20] can be of the order of the 
cavity linewidth. The idea is as follows: Consider the 
MSI operated at a point where the transmissivity r of 
the effective mirror is close to zero; r will then depend 
sensitively on the membrane displacement x. Combining 
this compound 'MSI-mirror' with a perfect mirror in its 
dark port will result in an effective FPI whose linewidth 
depends on x dominantly via t and not via the change 
in the (effective) cavity length L, see Fig. 1. In con- 
trast to a true FPI, therefore, g^ and g^ have a different 
functional dependence on x. This feature gives rise to 
a topology (i) where the optomechanical interaction can 
be tuned at will between being strongly dissipative or 
dispersive, and (ii) which can realize dissipative optome- 
chanics leading to ground-state cooling of the mechanical 
oscillator via quantum noise interference, as discussed in 
Ref. [24]. The distinctive signature of this interaction, 
an asymmetric Fano lineshape, is observable in the spec- 
trum of the cavity output field. We also report on the 
suitability of using this system for (iii) sensitive position 
transduction. 

This paper is structured as follows. We shall first de- 
scribe the physical model and write down the resulting 
Hamiltonian and input-output relation. These are then 
used to derive the equations of motion for the cavity field 
and oscillator motion. The resulting dynamics is solved 
to obtain the steady-state mechanical occupation number 
analytically in the weak-coupling limit, and also numer- 
ically in both weak- and strong-coupling regimes. 



I. MODEL 



dispersive-only Hamiltonian that is considered in most 
works on optomechanical systems. Before we proceed to 
make the connection of Hamiltonian (1) and the physics 
of a MSI [31, 32] we note that the standard cavity- input- 
output relation must be generalized to accommodate the 
effect of dissipative coupling (see Appendix): 
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This relation provides a boundary condition connecting 
light leaving the effective FPI with light entering it, and 
the intracavity dynamics. It can be identified with the 
familiar relation Cout — Oin = \/'^i^{i) taking k{x) = Kc + 
gf^x and truncating the square-root to first order in x. 

The generalized optomechanical Hamiltonian (1) is re- 
alized in the MSI shown in Fig. 1 operating close to a 
dark-port condition, i.e., when most light is directly re- 
flected at the first beam splitter (BS). The whole system 
can then be described as an effective Fabry-Perot inter- 
ferometer of length L, operating in the good-cavity limit, 
formed between the perfect end-mirror Mi and an effec- 
tive end-mirror M.\ 2L is the length of the Sagnac mode 
M1-BS-M2-M3-BS. The reflectivity p and transmissiv- 
ity r of A^ depend on the complex reflectivity R (r) and 
transmissivity T (i) of BS (of the micromirror), as well as 
the displacement (5 of M from the mid-point of M2-M3: 

p = -(E^re^**' + T^re-^'-J + 2i?rt)e-*"s' , and (3a) 

"■s*. (3b) 



T = [{RT*e''^ - c.c.)r - {\R\'' - \T\')i]i 



Assuming the close-to-dark-port condition \t\ <C 1, 
a quantization along the standard routes of cavity 
QED (see Appendix) yields a Hamiltonian of the form 
in Eq. (1) with the usual result Kc — —c/{2L) ln]p) « 
c]t]^/(4L), as well as dispersive and dissipative couplings 



The Hamiltonian including both dispersive and dissi- 
pative effects in an opto- or electromechanical [24] system 
can be written (see Appendix) 
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[a; ,a — H.c. 



This Hamiltonian is written in a frame rotating at the 
optical frequency of some driving fleld Wd! A = Wd — Wc 
is the detuning from cavity resonance (ojc) whose annihi- 
lation operator is a. x = {h-\-o^)/ \/2 is the dimensionless 
displacement, and h is the annihilation operator for the 
mechanical oscillator of frequency Wm. The second line 
in Eq. (1) provides the coupling of the cavity mode a to 
the modes a^ of the external field coupling into and out 
of the effective FPI, and will give rise to the finite width 
Kc of the resonance. The two terms proportional to x de- 
scribe the shifts of the cavity resonance and width with 
mechanical displacement, with strengths characterized 
by gu>.K, respectively. This Hamiltonian generalizes the 



gu, = {oo,Xo/L)[{\R\^ - JTI^) +TCOs(argt)] , (4a) 
g^ = i\T\e''^{uj,xolL) [2RT -f pcos(argt)] , (4b) 

where « close to resonance. These results guarantee 
that the values of g^j and (7k can be controlled indepen- 
dently by choosing S {i.e., positioning the membrane) and 
the reflectivity of the central beam splitter BS appropri- 
ately. The need for a sharp resonance demands \R\ ~ \T\. 
We note, however, that \R\ ^ \T\ is required to be able 
to set gcj = with \t\ > 0. In the following we will spe- 
cialize to this most interesting case of a purely dissipative 
coupling. We will provide experimental case studies be- 
low and show that strong coupling, where [aUgKl 3> Kc 
(a being the intracavity amplitude), can be achieved for 
a moderate driving power of a few hundred- |J.W. 

The next step in our investigation is to derive the 
Heisenberg equations of motion for a, x, and p — [b — 
W)/{i\/2^. The Hamiltonian (1) implies a full and rich 
nonlinear dynamics [33], including such effects as bista- 
bility and self-induced oscillations, but in this new con- 
text of dissipative optomechanics. Here we will focus on 
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FIG. 2. Calculated steady-state occupation number for sys- 
tem II. Throughout this plot, A = a;in/2 and Pin = 10 nW. 
Inset: Normalized backaction force noise density, illustrating 
its Fano profile, when LJm ~ 0.6kc and A = a;m/2; the noise 
density at to = —uum is zero. 



the linear dynamics by assuming a strong classical driv- 
ing field Oi^, which allows us to write the linear equations 
of motion for the fluctuations around the steady state 
as (see Appendix) 



a 



^iA- Kc)a- ^/2K^ain - gn {a-m/V'^Kc + a)x , (5) 



igK/\/2K^[(ain*a + aa-n) - H.c] 



(6) 



and X = uJraP, with the mechanical motion having damp- 
ing rate Km- a is the coherent part of the cavity field; 
and ^ = ^^, which models Brownian-motion-type noise 
acting on the mechanical oscillator, is assumed to obey 
{i{t)iit')) = (2nth + l)'5(t - t'), where nth is the ther- 
mal occupation number in the absence of driving, and 
(^) = 0. The condition for the c-number component OiTI" 
of the input field to be 'large enough' is obtained by eval- 
uating the contributions of both components of the input 
field to (a^a). This linearization condition can be stated 
as loi^P > [(|A|+CLim)^-|-K^]/(2Kc) (see Appendix), and 
will be assumed to be satisfied in the following. In the 
weak-coupling limit, when I^kPIoi^P <C 4k^, we assume 
that the cavity field follows the mechanical motion adia- 
batically and solve the linearized dynamics to obtain the 
steady-state mechanical occupation number: 
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(7) 



In this equation Wm and R,^ are the optically-shifted me- 
chanical oscillator frequency and damping rate, respec- 
tively, whose expressions are rather involved and will not 
be reproduced here (see Appendix). The preceding re- 
lation is valid even for moderately strong input powers, 
as can be shown by comparing this analytical result with 
that obtained by solving the above equations of motion 
exactly using the method in Ref. [34]. In the 'cryogenic 
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FIG. 3. (Color online) Full numerical solution for the occu- 
pation number (Vb) for (a) system I, and (b) system II. 



optomechanics' limit, where Kc 3> Km ^ 'i-th'^m and 
Wfn w ojnij this result simplifies to 



(6t6) 
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4Km A2 + k2 (A - Wm)^ + «:? 



(8) 



This expression for {b^b) deserves some comments. We 
recall that in dispersive optomechanics, cooling is op- 
timized when the upper motional sideband is strongly 
enhanced (A = — Wm) and the lower sideband strongly 
suppressed {cOm 3> Kc)- In the present case, Fano-like in- 
terference can be observed in the backaction force noise 
spectrum Spp,{uj), cf. the inset of Fig. 2, where F [the sec- 
ond line in Eq. (6)] is the force operator acting on the mi- 
cromirror motion. This resonance neutralizes the lower 
sideband when A — ujj^/2; the Fano lineshape also means 
that optimal enhancement of the upper sideband requires 
not Wni ^ Kc — indeed, that would be disingenuous — but 

opt _ 



W„ 



2V13-5, 



' 0.6kc, which is much less 
demanding than the sideband-resolved condition. Fig. 2 
shows how the mechanical occupation number, at the 
optimal detuning A = Win/2 and with an input power 
Pin = 10 nW, changes as the ratio uim/nc is varied, and 
is minimized when w,,, sa w°f ' . 



II. COOLING 

Let us now turn our attention towards predicting the 
cooling performance of the model investigated above. We 
shall call 'system I' the following: mechanical oscillator 
effective mass m — 100 ng, frequency w^ = 27r x 103 kHz, 
and quality factor Q = ui^n/i^Km) — 2 x 10^, such that 



the zero-point fluctuation is xq = ^Jhji^hrujj^ ~ 1 fm; 
|i?|2/|r|2 = 0.486/0.514; |r|V|t|2 = 0.362/0.638; driving 
wavelength Ac — 1064 nm; and L — 7.5 cm; which are ex- 
perimentally realizable [32] and yield Kc = 27r x 196 kHz, 
5«, « -27r X 0.1 Hz [g^jxQ « -27r x 79kHz/nm). We 
limit the input power to the regime where heating from 
the power absorbed is not the dominant process. An in- 
put power of 10 mW corresponds to an effective coupling 
strength G = |a||5K| ~ 0.1 Kc- By starting from an envi- 
ronment temperature Tgnv = 300 K we can decrease the 
occupation number by over three orders of magnitude, as 
illustrated in Fig. 3(a). 

Consider next a hypothetical, but still physically realiz- 
able, situation ('system H') where M has a smaller mass 
(m = 50 pg [29]), higher mechanical quality (Q = 1.1 x 
10^ [35]), and higher reflectivity (|r|V|t|2 = 0.818/0.182, 
possible by patterning the membrane [36]), where BS 
is more balanced {\R^ l\T\^ = 0.496/0.504, Kc = 27r x 
59 kHz), and thus g^, ~ — 27r x 2.6 Hz {qk/xq ~ — 27r x 
65 kHz/nni) . For this system, we can achieve the strong- 
coupling condition G > Kq- If we also assume cryogenic 
operation at Tenv = 0.3 K, we can see from Fig. 3(b) that 
ground-state cooling is possible, despite the poor reflec- 
tivity of the mechanical oscillator, and that Wm ~ Kc- 

Position transduction. — Optomechanical systems are 
one promising approach towards extremely sensitive po- 
sition transduction [1]. The coupling of a mechanical 
oscillator to multiple cavity modes was recently explored 
in Ref. [37]; a common feature of 'multiple light mode- 
single mechanical mode' systems is a Fano-like profile in 
Spp{uj) (Fig. 2 inset). The corresponding anti-resonance 
allows one to reach the SQL for measurement impreci- 
sion at a significantly lower input power than a 'single 
light mode-single mechanical mode' optomechanical sys- 
tem. Indeed, let us quantify the achievable resolution 
of a position measurement by the ratio Af/S, where the 
'noise' Af"^ is the contribution to the symmetrized homo- 
dyne output spectrum evaluated at a; = Wm, So^it^uj^^), 
due to Sin, and the 'signal' S^ is that due to ^, normalised 
to the free mechanical motion noise spectrum S'fi.ce('^m): 



5'out(Win) = 7V^ + 5^ S'frco(Wm) ■ 



(9) 



At A = and under identical conditions, it can be shown 
that J^/S reaches the same lower bound (the SQL) in 
both dissipative (at a power P^) and dispersive (Pai) 
cases, but with P„ = {2K/uJra)'^Pu- In the sideband- 
resolved regime, therefore, one can obtain significantly 
better position resolution at low powers in comparison 
with the dispersive case. 



Ref. [38], the membrane temperature rises by ca. 60 K 
at Pin = 1 mW. The resulting temperature rise has a 
significant effect on the base occupation number of the 
micromirror, but still allows strong cooling of the mi- 
cromirror motion. It is a feature of our topology that 
the ideal position of M corresponds to both where Pabs 
is greatly reduced, and where the competing dispersive 
optomechanics is switched off. 

Significantly, we note that this situation does not per- 
sist in the case of the 'membrane-in-the-middle' geom- 
etry. In a single-transverse-mode model for this latter 
situation, the dispersive and dissipative optomechanics 
cannot be independently turned off, and are both zero 
at the nodes of the cavity field; this leads to a stronger 
restriction arising from the power absorbed. Moreover, 
dispersive (dissipative) optomechanical cooling requires 
A < (A > 0); most treatments of this geometry do not 
include the dissipative optomechanics component of the 
dynamics [39], and (at least within a single-mode model) 
may therefore over-estimate the cooling efficiency in the 
non-sideband-resolved regime. 



IV. CONCLUSIONS 

In this paper we have presented an experimentally fea- 
sible realization of an optomechanical system that can be 
fully tuned between strong dissipative and dispersive dy- 
namics. The cooling mechanism we presented works best 
in the non-sideband-resolved regime, unlike the usual dis- 
persive case. For an existing set of experimental param- 
eters, we predict a strong cooling effect arising from a 
Fano-type resonance in the backaction force acting on 
the mechanical motion. In the case of more optimistic pa- 
rameters, we predict ground-state cooling of the mechan- 
ical motion. In the opposite, sideband-resolved, regime, 
our system promises significantly improved position mea- 
surement resolution. The usage of the proposed imple- 
mentation of a dissipative optomechanical coupling for 
nonlinear (quantum) dynamics [26, 27], and for pulsed 
dynamics [28] remain to be explored. 
Moreover, we believe that the topology of a MSI with 
a highly refiective mirror in its dark port — the signal 
recycling configuration used in the context of gravita- 
tional wave detectors [40, 41] — and the resulting dissipa- 
tive intracavity dynamics provides attractive possibilities 
for other fields in cavity QED, especially in combination 
with single atoms or ensembles of cold atoms. 



III. COMMENTS 

We set gu} — early on in this paper, which can be 
realized by placing M at a point where the field inten- 
sity surrounding it is close to minimal, greatly reduc- 
ing the power absorbed Pabs by the membrane. For the 
parameters in Fig. 3(a), using the 'thermal link' from 
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FIG. 4. (a) Schematic layout of the Michelson-Sagnac inter- 
ferometer. The z-axis is defined as having its origin at Mi. 
X is the operator representing the displacement of M from its 
equilibrium position, (b) The optics in the dotted box in (a) 
form an effective mirror, A4. 
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FIG. 5. (a) The four field modes interacting with the mi- 
cromirror. (b) The eight field modes interacting with the 
beamsplitter. Ci,2 and -Di,2 are identical, up to a phase shift, 
to Ci 2 and -Di,2, respectively, (c) Phase shifts incurred in 
each of the three sections of the interferometer, at the cavity 
frequency uic- 



Appendix A: Field inside the interferometer 



Z3. For a beam of cross-sectional area A, the frequency- 
dependent normalization factor is Af{uj) = 



hu 



With di^ being the annihilation operator for the mode 
at frequency w, we can thus write the field Hamiltonian 
ih=l) 



H ^ duj 



(A3) 



In this half-space, we can solve for all the field mode 
amplitudes, ^1,2, B2, Ci_2, C'l 2, D1.2, and D[2, in terms 
of Bi by, e.g., using the transfer matrix formalism [42, 
43] . Upon writing this solution, it becomes apparent that 
the optics inside the dotted box in Fig. 4 can be treated 
as a single effective mirror M, having reflectivity and 
transmissivity 
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R*Txe 



-iS 



\R\' 



\T\^)i] 



Ji'^-x) 



(A5) 



where i' is the phase-delay Mi-BS, 6 — 9 — (f> is the arm- 
length difference for the loop in the interferometer, cf. 
Fig. 5(c), and x — a'i'g(i). In this notation, we can now 
write, e.g., 



_ T e^p[i{a + X + I')] p 
1 — pexp[^(c^ + X + 1^)1 



(A6) 



for convenience we have defined a = 9 + cf). The coeffi- 
cient linking the "intracavity" field amplitude B2 to the 
one of the incoming field Bi is completely equivalent to 
the familiar spectral response function of a Fabry-Perot 
cavity, see e.g. Ref. [44]. 



Consider the Michelson-Sagnac interferometer topol- 
ogy shown in Fig. 4, where all the fixed mirrors are per- 
fect. The micromirror, assumed lossless, has amplitude 
reflectivity r and transmissivity t. The corresponding 
parameters for the lossless beamsplitter are R and T, re- 
spectively. The electric field for the half-space defined by 
z > can be written as 



E^{z)= / dujM{uj)ui^{z)au: , 
Jo 



(Al) 



where the mode function 



Uuj{z) 



Aze"**-'^ + B2e'''^ 


0< z <zi 


Cie-'^' + Die'''' 


Zl < Z < Z2 


D2e-'''' + C2e'''' 


Z2 < Z < Z3 


Bie-''''' + Aie'''^ 


Z3<Z 



(A2) 



with reference to Fig. 5. Starting from z = 0, BS is at 
z = Zl, M at z = Z2, then BS is encountered again at z = 



Appendix B: Interaction Hamiltonian 

In this section, we will derive an effective Hamiltonian 
Hint that describes the interaction between the motion of 
the micromirror and the electric field in the regime where 
the displacement x of the membrane from its average 
position (fixed by S) is small on the scale of a wave length, 
i.e. in the Lamb-Dicke regime. The Hamiltonian for the 
radiation-pressure interaction in linear order of x is [45] 



Hint ^ Aeoiv^ - l)[E- {z2 - 
-E-(z2 



)E^ 



(Z2 - 

'l)E^Z2 



)]x, (BI) 



with membrane thickness d and refractive index rj. We 
now make the simplification that the membrane is in- 
finitesimally thin, which corresponds to taking d — > 
and ry — > (X) such that r/^d is constant. This allows us 
to express Hint in terms of left (right) propagating fields. 



Field solution, Eq. (A6) 

Good-cavity limit 

Field Hamiltonian — 



Lamb— Dicke 
expansion 



Lamb-Dickc 
expansion 



Eq. (Bl) We can then write the interaction Hamiltonian (B3) as 



Good-cavity limit 

i 

--^Eq. (C3) 



FIG. 6. Schematic of the alternative paths to obtaining the 
interaction Hamiltonian. Here, we follow the clockwise path 
shown by the solid arrows, but an equivalent formulation fol- 
lows the dashed arrows. 



E- 



\r)^ through Wl(r)^„(2:): 



(B2a) 
(B2b) 
(B2c) 

(B2d) 



which are expressed in terms of Bi and B2 in an- 
ticipation of taking the good-cavity limit, and where 
Z2 = hmc;^o(-Z2 ± 5)- With this notation, we have (cf. 
Ref. [43]) 

iJint = 2Aeo [E^ {z2 )El{z2)+ E^ (z^^ ) E+ (z^- ) 

)E+{z+)-EAz+)E+{z+)]x. (B3) 



-e; 



Appendix C: Good-cavity limit 

We now go to the good-cavity limit in our system. This 
condition, which is often termed the 'dark port' (or 'dark 
fringe') condition [31], corresponds to having |p| w 1. 
Under such conditions, the system acts as an effective 
Fabry-Perot interferometer with one perfect end-mirror. 
Ml, and the effective end- mirror A4. In other words, 
|i?2|^ exhibits a Lorentzian behavior, with the peaks oc- 
curring at the poles of the denominator of Eq. ( A6) , at 

Wc = ^ [arg(p) -icT + x + i') + 2nn] , (CI) 

with n being an integer, and having a half-width at 
half-maximum (HWHM) n^ = -c/(2i) ln|p|, L = {a + 
i')c/{2u!c) being the length of the interferometer. 

In the following we will concentrate on one particu- 
lar cavity resonances, and consider n to be fixed. For 
this single mode cavity field we introduce the bosonic 
field annihilation operators a and reassign cli^ to be the 
annihilation operators for the bath that the cavity field 
couples to; these operators satisfy the usual commutation 
relations 



+ 9o\j^(p j duj ilia + \l.c\b + P), (C3) 

where for clarity we have defined 

a = 2 [(|i?|2 - |T|2) -f re-"" cos(x)] , and (C4) 

/? = -2 [2RT + pe-"" cos(x)] e*('^+^+'') , (C5) 

We may relate a and /3 to g^^.K from the main text by: 



9i^ = .9o(a/V2) , and .g^ = -igaPd 



4:KrL 



-Wo\t\(3. 

(C6) 

which are assumed to be slowly varying over the fre- 
quency range of interest. The couphng constant go, which 
is dimensionally a frequency, is defined as go — ujcXo/2L, 
and b represents the annihilation operator for the motion 
of the membrane. In terms of b and the zero-point fluctu- 
ations xq, we can write the displacement of the micronrir- 
ror from its equilibrium position as x = 2:0(6 + b^)/\/2. 
(Note that in the main text we used a dimensionless x.) 
In the above formulation, we first computed the Hamil- 
tonian to first order in x and then took the good-cavity 
limit. An equivalent formulation would be to first take 
the good-cavity limit, with the linewidth of the cavity de- 
pending parametrically on the position of the membrane, 
and then expanding this linewidth to first order in x (cf. 
Fig. 6). We also remark that the analysis can be per- 
formed without taking the good-cavity limit, resulting in 
equivalent results but less transparent expressions. 



Appendix D: Equations of motion 

We can express the full Hamiltonian H as a. sum of 
the free cavity He, free micromirror H^, free field Hf, 
field-cavity interaction 7Jf_c, and cavity-micromirror in- 
teraction Hint Hamiltonians: 



H — Hc+ Hj-a + Hh+ Hh-c + Hint 

Hf. — LUcO, a , 

Hni =UJmb^b, 



Ht = i dij LO al.a,. 



Htr = i 



doj [aia - a^aS) 



(Dla) 

(Dlb) 
(Die) 

(Did) 
(Die) 



[a(.j, aL] = (5(0; - J) , and [a, a^] = 1 



(C2) and Hint is given in Eq. (C3). 



1. Cavity bath: Input— output relation 

We will first calculate the equation of motion for a^ : 



: *^\/i^(^+^') 



a. (D2) 



This equation can be formally solved (see Ref. [46, §5.3]) 
and, in the weak-coupling approximation, integrated over 
frequency to give 



2 ^i5o/3y£(&+fo^) a, (D3) 



where we have defined the input field 

am = ^ / dw e'^'^^'-'^^a^^o , (D4) 

V Stt J 

to being some remote initial time at which d^^ = d^^.o- 
We can similarly define an output field at some distant 
future time ti: 



^ 


1 / 


dwe-'"(*-*i)<i. 


flout 


V2^J 


and infer 






1 f 


= flout — 


/f-..o/3/|(^ + St) 


— / do; a^ 

V^ J 



(D5) 



(D6) 

Finally, Eqs. (D3) and (D6) give us the input-output 
relation 



flout flin 



2c 



a. (D7) 



The input field operators model a heat bath at zero tem- 
perature, and therefore satisfy the relations 

{aUt)al{t'))=Sit-t'), and (D8a) 



(flLWain(i'))=0. 



2. Cavity field 



(D8b) 



The equation of motion for the cavity field, correct to 
linear order in terms of the form {b + b^), reads 

a = —i[cJc + goOi{b + fo' )]d 



«e + 5oIm{/3}\/^(6 + 6^) 



" 2c 



in. (D9) 



Let us now define the effective cavity resonance frequency 
and linewidth, 

ijc^ujc. + goa{b + P), (DlOa) 

k,^K,+golm{l3}^'^{b + b^), (DlOb) 

respectively. Algebraic operations on Re are to be un- 
derstood in terms of series expansions in powers of e = 
y^KcL/c, correct only to first order in this parameter, 
e.g., 



V^ = v^ + igo Im{/3} J-ib + b'f). 



(Dll) 



We can simplify our expressions considerably if we intro- 
duce an effective coupling parameter F, which is dimen- 
sionally a frequency, through the relation 



(D12) 



2r^goRe{(3}\jg{b + P). 



where algebraic expressions involving F are to be in- 
terpreted similarly to those involving R^', in particular, 
2F = (\/2r) = to first order in e. In terms of these 
quantities, we can rewrite the input -output relation as 



flout fli] 



{V^Rc -iV2f)d 



(D13) 



Moreover, the equation of motion for d takes the simpli- 
fied form 



-{iuJc + Kc)a — (v 2kc + i V 2F)c 



Micromirror motion 



(D14) 



The motion of the micromirror is modeled as a single- 
frequency harmonic oscillator. The equation of motion 
for b reads 



b= -ito^b-igQaa^a- igo\/ :^ (/^dj^d -f /3*d'''din) 



lgo^|^{Pa + P*a^)■ (D15) 



We notice immediately that the last three terms in 
Eq. (D15) do not appear in the equation of motion for x] 
they describe a Hermitian force operator 

F = —V2mujniXogoa a' a 

- muJmXogoJ — {l3al^a + I3*a^ain) 



.xogo\J^{Pa + f3*a^). (D16) 



Defining the momentum quadrature of the micromirror 
motion as p = muj„iXo{b — 6')/(v2i), and introducing 



the motion-heat-bath couphng through the Brownian- 
motion Langevin operator ^ we obtain x = —p/m, and 



P 



-mujl^x ~ 2Knip + 2y/K^^mujiTiXo^ + F , (D17) 



where (|) = and (|(f)l(^')> = (2nth + l)S{t - t'), with 



"•th 



o?iW„,/(feBT,„v) _ 1 



(D18) 



being the thermal occupation number, in the absence of 
driving, at an environment temperature Tcnv 



Appendix E: Linearized dynamics 

The equations of motion derived in the previous sec- 
tion must be hnearized with respect to both a and ajn. 
In the case of the former, the criterion for hnearization 
can be stated simply: |ap ^ {a^a) = 1; in the case of Sin 
the situation is somewhat more involved. To see what 
the relevant criterion is let us examine the relevant con- 
tributions of iiinX and Oi^x to {a^a). Consider first, in 
the weak-coupling approximation, 

k{t) = (zA-A.e)aW+7«mW(^oe-^"'"*+&J,e^"»*) , (El) 

where we extract the harmonic time dependence, at fre- 
quency Wni, from 6, and where 7 is some constant. Let 
us again assume that 60 varies slowly. We can Fourier- 
transform this equation, compute (a^[a;]a[a;']), and re- 
verse the Fourier transforms to obtain 



We will now look at the c-number contribution. 



{aHt)a{t))^^-^{{blbo) + {bobl))- (E2) 



h = (iA - Kc)a + ja^Jboe'"'^"^* + SJe'""'*) . (E3) 

This equation can be solved, for the assumed case of weak 
coupling and slowly varying a^ and 69 • Upon computing 
(d^(i)a(i)), we can state the linearization criterion as 



.|2 ^ i\A\ + u;^y + ^ij 



(E4) 



Let us now calculate (b^b) in the steady state. We start 
with the two coupled linearized equations of motion for 



a = and imaginary /3: 
a — (jA — Kc)a — \/2Kcain 



/2^a)(6-f-&^), (E5a) 



and 



^{iuJn, + Kn,)b- \/2Km b;, 



igoP 



Oin a-VaalJ - H.c 



(E5b) 



In the adiabatic and small-coupling limit, Eq. (E5a) is 
formally solved by assuming that b varies slowly on the 
relevant timescale. The solution to this equation is then 
substituted into Eq. (E5b) to yield a differential equation 

relating b to the input fields and, in the rotating-wave ap- 
proximation, exclusively to b. This equation of motion for 
b exhibits optically-induced damping of the micromirror 
and the optical spring effect; the micromirror oscillation 
frequency is Wm -I- uJq , where 



-5o^l/3p|l^ 



1 



A2 



2AKg + (A + ^,n)(<-A2) 
(A + a;,n)2 + k2 

2A^g-KA-c^„,)(A.g_^2) 

(A - uj^Y + kI 

its damping constant is K,n + k^, with 



(E6) 



5ol/3|'||aI. 



Kjc 



(2A + uj- 



A2 -)-k2 
2 



(A-Kw„ 



Kt. 



{2A ~Lu^f 
(A - LU^^Y + ^l 



(E7) 



After some algebra we can finally show that 



2c A2 + k'I Km + Ko 

X [(A - CJm - CJo)^ + {Hc + Km + Ko)^] 
X {(Km + Ko) [A2 + kI + Kc{Kra + Kq)] 

+ Kc(2A - W,n -^o)"^} 
- "-thKm/ (Km + Kq) . (E8) 
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